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ABSTRACT 

We present high resolution (R~170,000) Kitt Peak National Observatory Coude 
Feed telescope observations of the interstellar K I 7698A line towards 5 multiple star 
systems with saturated Na I components. We compare the K I absorption line profiles 
in each of the two (or three) lines of sight in these systems, and find significant dif- 
ferences between the sight-lines in 3 out of the 5 cases. We infer that the small scale 
structure traced by previous Na I observations is also present in at least some of the 
components with saturated Na I absorption lines, and thus the small scale structures 
traced by the neutral species are occurring at some level in clouds of all column den- 
sities. We discuss the implications of that conclusion and a potential explanation by 
density inhomogeneities. 
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1. Introduction 

There is convincing evidence from recent radio 
and optical observations that the diffuse inter- 
stehar medium (ISM) exhibits significant subparsec- 
scale variations down to limits of a few AU (Frail 
et al. 1994; Davis et al 1996; Meyer k Blades 
1996; Watson & Meyer 1996). On somewhat 
larger scales, strong variations in interstellar Na I 
profiles are observed towards multiple late-type 
stars in globular clusters (Langer, Prosser, & 
Sneden 1990; Bates et al. 1995). In some cases, 
the inferred densities of these structures far ex- 
ceed the nominal diffuse cloud values and ap- 
proach those of molecular cloud cores. For exam- 
ple. Frail et al. (1994) found pervasive variations 
in the H I opacity on scales of 5-100 AU utiliz- 
ing multiepoch observations of 21-cm absorption 
toward high velocity pulsars, which implied den- 
sities of uh 10^-10^ cm~^. It is not clear how 
such small, dense structures can arise or be main- 
tained in the low pressure environments of diffuse 
clouds. 

One difficulty with using the Na I D lines as 
a tracer of small scale structure is the increas- 
ing saturation of these lines as one goes to larger 
H I column densities. Thus the bulk of the gas 
in heavily reddened sightlines cannot be sampled 
using the D lines, and one is generally limited 
to identifying variations in the weaker wings of 
the line. A potential probe of these sightlines 
is the A7698A interstellar K I line. The lower 
cosmic abundance of K coupled with differences 
in the photoionization cross section leads to K I 
columns that are typically a factor of ~50 be- 
low those of Na I in diffuse interstellar clouds 
(Hobbs 1974). Thus one can use K I to see if 
the small scale structure phenomenon extends 
from the very diffuse components that have been 
well sampled using the Na I D lines (Watson & 
Meyer 1996; Meyer & Blades 1996) to highly red- 
dened sightlines such as HD 206267 (this paper, 
E{B -V) = 0.50). 

In this Letter, we present Kitt Peak National 



Observatory (KPNO) Coude Feed echelle obser- 
vations of the interstellar K I absorption towards 
members of 5 multiple star systems which have 
saturated interstellar Na I profiles. These ob- 
servations allow us to trace the small scale vari- 
ations previously observed using Na I to much 
more heavily reddened lines of sight. In addi- 
tion, we present a simple method for estimating 
the possible differences in density in the clouds 
using the observed K I column variations. 

2. Observations 

The observations were obtained in June 1998 
with the 0.9m Coude Feed telescope and spec- 
trograph at Kitt Peak National Observatory us- 
ing camera 6 in echelle mode and a Ford 3000 
X 1000 pixel CCD chip. The resolution of the 
data was measured using the ThAr lamp emis- 
sion lines and is ~ 1.75 kms^^ at the location of 
the K I A7698.974A hue. A total of 5 systems 
were observed: 4 binaries or common proper mo- 
tion doubles and one triple system (see Table 1). 
Observations of the stars a Aql and a Leo were 
obtained template for dividing out telluric 
absorption in the vicinity of the K I 7698A line. 
For the majority of the stars, exposures were 
taken on different nights at different grating tilts 
to reduce the effect of any flaws in the CCD chip. 
However, due to poor weather conditions obser- 
vations of p Oph were limited to a single night 
and only a single exposure was obtained for HD 
206267 D. 

Reduction of the data was done with the 
NOAO IRAF0 echelle data reduction package. 
The individual frames were first bias-subtracted 
and flat-fielded, then the scattered light was re- 
moved. For all stars except p Oph the individual 
orders were then extracted and cleaned of cos- 
mic ray hits. Since both members of the binary 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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p Oph were in the slit, only the portions of the 
combined stellar profile that were dominated by 
either one of the two stars were extracted. All of 
the resulting one-dimensional spectra were then 
wavelength calibrated, shifted to heliocentric co- 
ordinates, summed, and finally continuum fitted 
using low order polynomials. Figures 1 and 2 
show the final K I line profiles. The resulting 
S/N-ratios for these observations varies widely, 
from as low as ~20 for the faint star HD 206267 D 
to as high as ~ 130 for fP Sco. Column densities, 
b- values, and relative velocities of the various 
components were derived by profile fitting with 
the programs xvoigt (Mar &: Bailey 1995) and 
f its6p (Welty, Hobbs, & York 1991). The wave- 
length and oscillator strength of the K I 7698A 
line were taken from Morton (1991); the inclusion 
of hyperfine splitting (Welty, Hobbs, & Kulkarni 
1994) in the fits does not significantly alter the 
derived columns or b- values. The column den- 
sity, b-value, and velocity of the various com- 
ponents were free parameters in the initial fits, 
while in the final fits the velocities of the com- 
ponents were fixed since the velocity differences 
between members of a system were smaller than 
the uncertainties. The resulting column densi- 
ties are listed in Table 2; the listed uncertainties 
include contributions due to signal fluctuations, 
continuum placement errors, and the errors in- 
duced by uncertainty in the b- values. It must 
be noted that if there is unresolved component 
structure that the column density may be higher 
than that listed. In particular, higher resolution 
observations of p Oph A suggest that there are at 
least three distinct components making up what 
is fit here as a single main component along this 
line of sight (Welty 1998). 

3. Discussion 

As seen in Figures 1 and 2, significant differ- 
ences in the absorption profile were seen between 
the interstellar K I lines towards the stars in 3 of 
the 5 systems. The detection of significant varia- 
tion in the profiles of the K I within the saturated 



cores of the (previously observed) Na I profiles 
shows that the small scale structures occur even 
in regions of high column density. It should be 
noted that the detection of widespread variation 
in Na I and K I in multiple components make a 
circumstellar origin for the observed small scale 
structures significantly less tenable. For exam- 
ple, towards /3 Sco variations in Na I are ob- 
served at V ~ —24 kms~^ and in K I in the 
main component at v ~ — 10 kms^-*^, with the 
enhancements in Na I and K I columns occurring 
towards different members of the binary. One 
also sees variations in multiple K I components 
towards the HD 206267 system, but one must 
note that the HD 206267 system bears a strik- 
ing resemblance to the Trapezium (Abt 1986), 
and like the Trapezium is in a region of recent, 
active star formation. Thus much of the vari- 
ation towards this system may be in gas which 
is located relatively near the stars (although not 
circumstellar), and may be the result of interac- 
tions between the surrounding medium and stel- 
lar winds and/or associated with the H II region 
IC 1396 centered on HD 206267A (O'DeU et al. 
1993). While large differences were not detected 
towards p Oph (HD 147933/4) or HD 161270/89 
it is possible that significant variation could be 
hidden in the strong (possibly saturated) cores of 
these lines. In any case, the observed K I varia- 
tions suggest that the small scale structures pre- 
viously identified using Na 1 are truly ubiquitous, 
and occur even in highly reddened sightlines. 

We can use the observed variations in the K I 
to infer the corresponding Na I variations for 
comparison with previous optical studies of small 
scale structure. Hobbs (1974) identified a rela- 
tionship between the total line of sight column 
densities of K I and Na I: 

N(K I)~ 0.012 N(Na I)^-^, 

where N(K I) and N(Na I) are in units of lO^^cm"^ 
Assuming the above "typical" relationship be- 



3 



tween N(K I) and N(Na I), we see that the ob- 
served K I variations correspond to Na I varia- 
tions as large as ~ 1.3 x lO^^cm"^. Thus small 
scale variations occur over 3 orders of magnitude 
in Na I column density (Meyer & Blades 1996; 
Watson & Meyer 1996; this paper), and therefore 
represent more than a population of small, low 
column clouds. Furthermore, if we assume the 
"typical" relationship between N(K I) and N(H I) 
from Hobbs (1974), we can estimate the density 
nH(cm~^) in the structures responsible for the 
observed profile variations using the binary sep- 
aration by assuming that variations in K I trace 
variations in H I. Then the observed N(K I) dif- 
ferences imply Hff > lO^cm"^, similar to the 
densities inferred for a number of components 
in previous Na I studies (Meyer & Blades 1996; 
Watson &: Meyer 1996). Such densities are also 
only slightly less than the densities {nn ~ 10^- 
lO^cm"^) inferred by Frail et al. (1994) in their 
study of 21-cm absorption towards high velocity 
pulsars. 

The lack of observable changes in the Zn II col- 
umn density in components where the Na I col- 
umn varies towards the binary ji Cru suggest that 
at least some fraction of the variation detected in 
the Na I studies is due to density, temperature, 
or ionization fraction fluctuations, and thus not 
indicative of variations in the hydrogen column 
density between the lines of sight (Lauroesch et 
al. 1998). If we assume that the observed K I 
variations are due to similar fluctuations, we can 
use the measured column density differences to 
estimate (albeit somewhat crudely) the difference 
in density and temperature in these structures. 
The recombination rate for potassium is propor- 
tional to ~ Ue X T~°-^ (Pequignot & Aldrovandi 
1986); if we assume a neutral ideal gas in pres- 
sure equilibrium then the recombination rate will 
go roughly as ~ Ug x n^^. Typically it is as- 
sumed that the Ue/nn ratio is roughly constant 
in neutral interstellar clouds, with the dominant 
source of electrons being carbon atoms. Under 
this assumption we can then estimate required 



density contrast from the observed columns since 
N(K l)oc Ujj'^. We must note that as the den- 
sity increases and/or the temperature decreases 
in these clouds there will be a corresponding in- 
crease in the N(C l)/N(C II) ratio, and thus a 
corresponding decrease in the Uc/uh ratio. How- 
ever, if N(C I) « N(C II) in both clouds then 
the Ue/uH ratio will be roughly constant. In 
any case, we will initially assume Ue/nn is con- 
stant, and then identify (if necessary) any cases 
in which a correction is required. 

Table 2 lists the column densities from our 
profile fitting analysis of these lines of sight. 
Based on the measured differences in the column 
densities for the various components we have es- 
timated the density contrast Juh assuming an 
ideal gas with a constant Ue/nn, then: 

<5nH = [N(K l)stron9M'K I)^eafe]°-' 

where we have defined the strong and weak com- 
ponents such that N(K l)strong >N(K l)yjeak- 
Thus all values of (^njj will be > 1. Looking at 
Table 2, we note that the estimated density dif- 
ferences are generally quite small, generally less 
than a factor of 2. This shows how relatively 
small density and temperature fluctuations can 
give rise to large differences in the column densi- 
ties of neutral species. 

There have been a number of simplifying as- 
sumptions made in the preceding argument which 
may not be correct. First, we have assumed that 
the observed column density differences reflect 
density and temperature fluctuations and not H I 
column variations. Despite the lack of Zn II col- 
umn density variations along the lightly reddened 
// Crucis line of sight (Lauroesch et al. 1998), the 
radio observations of pervasive variations in the 
21 cm opacity towards pulsars and extra-galactic 
radio sources imply that some (or even all) of 
these structures may be associated with H I col- 
umn density variations. In addition, these struc- 
tures may not be in pressure equilibrium with 
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the surrounding medium - since we have no in- 
formation about the hfetimes of these struetures 
we cannot say whether they are stable structures 
or not. We have also assumed that the thermal 
gas pressure dominates in these structures, im- 
plicitly ignoring any magnetic field effects. Thus 
we may have underestimated the magnitude of 
the density fluctuations in the sightlines. On the 
other hand, we have assumed that the Uc/uh ra- 
tio is constant, another possibility is that there 
are pervasive fluctuations in the electron density 
on small scales or that the ratio decreases in the 
denser clouds due to the recombination of car- 
bon. Only by directly measuring the densities 
(ue and uh) in individual components towards 
multiple systems through the use of diagnostics 
such as the C I and C II fine structure lines can 
the nature of these structures be fully explored. 
Such ultraviolet observations will also enable the 
identification of any sightlines with H I column 
density enhancements similar to those inferred 
from radio studies. However, it must be noted 
that optical studies remain an efficient probe of 
small scale structures, and can be used to trace 
such structures on a variety of scales. 

It is a pleasure to acknowledge the support of 
the staff of KPNO, especially Daryl Willmarth. 
We would also like to thank Ed Jenkins, Eric 
Sandquist, Dan Welty, Don York, and the referee 
for their valuable comments and suggestions. 
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Fig. 1. — Observed interstellar K I profiles for 

both members of the binary star systems HD 
144217/8 ip Sco), HD 145502/1 {u Sco), HD 
147933/4 (p Oph), and HD 161270/89 taken us- 
ing the KPNO Coude Feed Telescope. Note that 
the vertical scale varies between plots. There 
are significant differences in the (relatively weak) 
profiles towards v Sco and /3 Sco, and a lack of 
apparent variation in the strong (possibly satu- 
rated) lines towards p Oph and HD 161270/89. 
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Fig. 2. — Comparison of the Na I D2 line profiles 
and K I profiles observed towards the HD 206267 
multiple star system. Note in particular the large 
number of variable components in the K I profile 
that are hidden within the broad saturated Na I 
profile. This shows the importance of K I as a 
probe of small-scale structure in regions of high 
column density. 
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Table 1: Stellar Data 



System 


Alternate 


ya 


E(B-V)^ 


Separation'^ 


Separation'^ 




ID 


(mag) 


(mag) 


(") 


(AU) 


HD 144217/8 


(5 Sco 


2.6/4.9 


0.21 


13.6 


2720 


HD 145502/1 


u Sco 


4.0/6.3 


0.25 


41.1 


7200 


HD 147933/4 


p Oph 


5.0/5.9 


0.48 


3.2 


740 


HD 161270/89 




6.2/6.6 


0.07 


20.6 


2880 


HD 206267 A/C/D 




5.6/8.4/8.0 


0.52 


d 


d 



"Yale Bright Star Catalog (HofHeit & Jaschck 1982). 

'Derived using the observed (B-V) and spectral type from the Yale Bright Star Catalog (Hoffleit & Jaschek 1982), and the 
intrinsic colors from Mihalas & Binney (1981) 

'^Projected separation between primary and sccondary(s) based on the spectroscopic parallax using the spectral types given 
in the Yale Bright Star Catalog (Hoffleit & Jaschek 1982). 

''Triple system with separations of A-C 11.7", A-D 19.9", and C-D 31.6". These angular separations correspond to projected 
separations of 14,300 AU, 24,300 AU, and 38,550 AU respectively. 
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Table 2: Column Densities 



System 


Velocity 


N(K I)^ 


b-value*^ 






(kms~^) 


(xlO^o cm-2) 


(kms~^) 


(cm~^) 


HD 144217/8 


-9.9 


12.3±1.0/16.6±2.4 


0.4±0.1 


1.2 




-11.7 


3.3±0.6/2.3±0.5 


0.4±0.1 


1.2 


HD 145502/1 


-9.3 


11.5±0.8/13.1±1.3 


0.7±0.1 


1.1 




-11.1 


7.1±1.0/2.6±1.1 


0.3±0.1 


1.8 




-13.3 


6.1±1.0/4.4±2.3 


(0.2) 


1.2 


HD 147933/4 


-7.8 


117±18/118±22 


1.3±0.1 


<1.3 




-12.5 


1.7±0.5/<1.3 


0.9±0.3 


>1.0 


HD 161270/89 


-16.5 


32.4±4.0/29.8±3.5 


0.85±0.15 


<1.2 


HD 206267 A/C/D 


-18.9 


53.5±4.7/86.8±10.6/55.3±12.6 


1.2±0.1 


1.3 




-17.2 


45.5±6.5/<3.6/54.3±14.7 


0.9±0.1 


>4.2 




-16.0 


20.5±5.5/20.5±8.2/<8.2 


(0.3) 


<1.6 




-14.5 


64.8±10.7/47.1±12.6/64.8±36.4 


0.3±0.1 


1.2 




-13.0 


28.2±2.7/36.0±7.6/48.5±12.5 


0.55±0.15 


1.2 




-10.8 


45.2±3.7/66.7±2.4/83.1±18.2 


0.6±0.1 


1.3 




-8.9 


6.9±1.7/4.2±2.6/6.9±3.5 


(0.2) 


<2.7 




-4.8 


<1.3/6.9±3.6/<2.0 


0.2±0.1 


>1.7 



"Column densities with Icr errors (or 2a upper limits) which include the uncertainties due to signal fluctuations, continuum 
placement errors, and the errors induced by uncertainty in the b-values. Values for the various members of the systems are 
separated by /'s. 

''Derived b-values, values in ()'s were fixed in the analysis and are less well determined. 

'^Estimated density contrast (see text). Note that for the HD 206267 system we give Snn only for HD 206267 A and C due 
to the relatively large uncertainties in the fits to the lower S/N-ratio data for HD 206267 D. 
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